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Nonlinear  two-dimensional  sail  theory 

Jean-Marc  Vanden-Broeck*1 

Department  of  Mathematics.  Stanford  University,  Stanford,  California  94305 
(Received  23  March  1981;  accepted  5  January  1982) 

Steady  two-dimensional  flow  past  a  sail  is  considered.  The  sail  is  assumed  to  be  supported  by  two  masts.  The 
flow  and  the  shape  of  the  sail  are  determined  as  functions  of  the  direction  a  of  the  flow  at  infinity  and  the 
Weber  number  A .  The  fuff  nonlinear  problem  is  formulated  as  an  integrodifferemiaJ  equation  for  the  shape  of 
the  sail.  This  equation  is  discretized  and  solved  numerically  by  Newton's  method.  Sail  profiles,  the  slack  in 
the  sail,  and  the  lift  coefficient  are  presented  for  various  values  of  a  and  A . 


I.  INTRODUCTION 

We  consider  the  deformation  of  a  two-dimensional 
sail  due  to  the  steady  potential  flow  of  an  incompres¬ 
sible  inviscid  fluid  around  it  (see  Fig.  1).  The  sail  is 
supported  by  two  masts  and  is  characterized  by  its 
constant  tension  o  and  by  the  distance  c  between  the 
two  masts.  The  fluid  has  density  p  and  velocity  U  at 
infinity.  As  we  shall  see,  the  shape  of  the  sail  is  de¬ 
termined  by  the  direction  a  of  the  flow  at  infinity  and 
by  the  Weber  number 

X  =  2  pcU2/0.  (1) 

This  problem  was  first  considered  by  Thwaites1  and 
Nielsen.2  These  authors  obtained  approximate  solu¬ 
tions  for  a  small  by  using  thin  aerofoil  theory.  Their 
work  was  further  generalized  by  Tuck  and  Haselgrove3 
and  Dugan.4  More  recently,  Vanden-Broeck  and  Kel¬ 
ler5  obtained  an  asymptotic  solution  for  x  small  for  ar¬ 
bitrary  values  of  a. 

In  the  present  paper  we  solve  the  fully  nonlinear 
problem  numerically.  In  Sec.  II  we  formulate  the  prob¬ 
lem  as  an  integrodifferential  equation  for  the  unknown 
shape  of  the  sail.  In  Sec.  ni  we  present  a  method  to 
solve  this  equation  numerically.  The  method  involves 
discretization,  which  converts  the  equation  into  a  set 
of  nonlinear  algebraic  equations.  Then  it  employed 
Newton's  method  to  solve  these  equations.  The  results 
obtained  are  discussed  inSec.  IV. 

Many  different  families  of  solutions  exist.  Tuck  and 
Haselgrove ’s’  stability  analysis  shows  that  only  one 
family  of  solutions  is  stable  for  a  small.  Therefore, 
we  present  only  numerical  results  for  this  family.  It 
is  found  that  for  each  value  of  a  there  is  a  maximum 
value  of  the  Weber  number  x  above  which  this  family 
fails  to  exist. 


II.  FORMULATION 

Let  us  consider  the  steady  two-dimensional  flow  of 
an  inviscid  incompressible  fluid  past  an  inextensible 
and  flexible  sail.  The  flow  configuration,  the  sail,  and 
the  coordinates  are  sketched  in  Fig.  1. 

The  leading  edge  or  luff  of  the  sail  is  attached  to  a 
mast  at  *  =  0,y  =0  and  the  trailing  edge  or  leach  is 
attached  to  another  mast  at  x  =c  ,y  =  0.  At  infinity  we 
require  that  the  flow  have  speed  U  and  direction  a. 

We  introduce  dimensionless  variables  by  taking  c  as 
the  unit  length  and  U  as  the  unit  velocity.  Let  the  shape 
of  the  sail  be  described  in  dimensionless  variables  by 
y  =  F{x ) , 0 «  is  1.  The  conditions  of  attachment  to  the 
masts  imply 

F(0)  =  F(1)  =  0.  (2) 

The  difference  between  the  pressures  on  the  two 
sides  of  the  sail  is  balanced  by  k(x)o,  where  k  is  the 
curvature  of  the  sail  and  o  is  the  constant  tension  in  it. 
By  using  Bernoulli’s  equation  we  obtain,  in  dimension¬ 
less  variables 

k(x)  =  (\/4  )(q‘(x-)l-?*(xr)2|.  (3) 

Here,  q'(x)  and  q‘(x)  are  the  flow  velocities  above  and 
below  the  sail,  respectively.  The  parameter  x  is  the 
Weber  number  defined  by  (1). 

It  is  convenient  to  introduce  the  complex  perturba¬ 
tion  velocity 

a  -ib  =  u  -  cos  a  -i(v  -sin  a) .  (4) 

Here,  u  and  r  denote,  respectively,  the  chordwise  and 
normal  components  of  the  velocity.  The  function  a  -ib 
is  an  analytic  function  of  z  -  x  +  iy  which  vanishes  at 
infinity.  Therefore,  we  can  write,  by  using  the  Plem- 
elj  formula 


The  integral  in  (5)  is  of  the  Cauchy  principal  value 
form. 

The  velocity  must  be  tangent  to  the  sail  on  both  sides. 
Thus,  we  have 


»’ Present  address:  Department  of  Mathematics  and  the  Mathe¬ 
matics  Research  Center,  University  of  Wisconsin,  Madison, 
Wisconsin  53706. 


(5) 

l - — - - - - - 

v\x)  =  u\x)F’{x),  0«r«l.  (9) 

Using  the  definition  (4),  we  can  rewrite  (8)  in  the  more 
convenient  form 

&•(*)  -  6-(x)  =  [«•(*)  - ,  (7) 

&•(*)  +  &’(*)  =~2  sin  a  4  [<i*(x)  4  a’(x)]F'(x) 

4  2  cos  aF’(x).  (8) 
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We  now  rewrite  (3)  in  the  form 

[l  +  £$]»»  +  4  [a*(X) -a'(x)Hl  +^'(*),]I2cob« 

+  a'(x)+a\x)]=  0.  (11) 

Finally,  we  impose  the  Kutta  condition  at  x  =  l  by  re¬ 
quiring 

«*(!)=<*  "(1)  •  (12) 


fig.  1.  Computed  profile  of  a  sail  with  a *=  t/6,  x  =  o.  82,  and  The  condition  (12)  is  clearly  not  satisfactory  when 
3  -  3. 0.  The  unit  length  in  y  is  c. 

F'(x)~(l  -xf  with  v  as  *  -  1 . 


By  substituting  (7)  into  (S)  and  taking  real  and  ima¬ 
ginary  parts  we  obtain 


&♦(*)  +  &*(*) 


n*(x)  +  a'(x) 


1  ('  la-(xl)-g-(xlm  +  F'(x,)»1IF(x,) -F(x) ) 
*  4  (x.-x^+lFtx,)-^*))2 


dx  j . 
(10) 


J 


However,  such  exceptional  cases  did  not  occur  in  the 
numerical  solution.  We  shall  consider  only  values  of 
a  between  0  and  n/2  so  that  the  Kutta  condition  is  im¬ 
posed  at  the  trailing  edge. 

For  given  values  of  x  and  a  relations  (8)-(12)  define 
a  nonlinear  system  of  integrodifferential  equations  for 
the  unknowns  F(x),  a*(*)>  a*(x)>  b’(x ),  and  6*(x). 

By  substituting  (7),  (9),  and  (10)  into  (8)  and  (11)  we 
obtain  the  following  reduced  system  of  equations  for 
the  two  unknowns  a*(x)  -a*(x)  and  F(x): 


•'0 


)-rwii 


(13) 


This  system  is  solved  numerically  in  the  next  section. 

From  the  solution  we  can  then  compute  the  slack  l 
in  the  sail  by  the  formula 

1=  /'{l  +  [F'(*)f},/Jrfx-l.  (15) 

■'0 

Following  Thwaites1  we  define  the  parameter  0  by  the 
relation 

^  =  sin2a/2/,/, .  (16) 

The  lift  coefficient  is  obtained  from  the  solution  by 
substituting  (7)  into  the  formula  (15)  given  by  Vanden- 
Broeck  and  Keller.9  Thus ,  we  obtain 

Ct=|2  [<T(x)  -a*(x)){l  +  [F'(x)]*]jfx  |  .  (17) 

III.  NUMERICAL  ANALYSIS 

In  order  to  solve  the  system  of  equations  defined  by 
(3),  (12),  (13),  and  (14)  we  introduce  the  mesh  points 

x,  =  (l-l)/(N-\),  1  =  1,. ..,N,  (18) 

and  the  corresponding  unknowns 

F,=F(x,),  1=1 . at,  (19) 

AJ  =a*(x,) -a~(x,) ,  1  =  2, . . .  ,N  .  (20) 

Relations  (2)  and  (12)  imply  F,  =  F„  =  0  andA^=0,  so 
that  there  are  only  2N  -  4  unknowns  F,  and  A", . 

We  shall  satisfy  Eqs.  (13)  and  (14)  at  the  N  -3  inter¬ 


mediate  mesh  points 

i /j  =  (*,  +  *,. i  )/2,  1—2 . N- 2.  (21) 

First,  we  compute  F'(x,),  F'fx,.,,,).  F"(x,.Wj),  F(x„ux), 
and  a*(x,.1/2)  -a'(x,.,/2)  in  terms  of  F,  and  A~,  by  four 
point  formulas. 

Next,  we  evaluate  the  integrals  in  (IS)  and  (14)  by 
using  the  trapezoidal  rule  with  the  mesh  points  xr 
The  mesh  points  x,  are  locally  symmetric  about  x,.,/} 
and  the  quadrature  formula  is  also  symmetric. 
Therefore,  the  contribution  from  the  neighborhood  of 
the  singularity  cancels  out,  permitting  us  to  evaluate 
the  Cauchy  principal  value  as  if  it  were  an  ordinary 
integral. 

In  doing  so  it  is  Important  to  notice  that  the  inte¬ 
grands  in  (13)  and(14)  behave  likex*,/,  as  x  —  0.  We 
handle  these  singularities  by  subtracting  the  singular 
part  before  integrating  numerically.  We  then  integrate 
the  singular  part  analytically. 

By  discretizing  (13)  and  (14)  we  obtain  2N  -6  nonlin¬ 
ear  algebraic  equations  for  the  2N  -  4  unknowns  F,  and 
A‘r  An  extra  equation  is  obtained  by  imposing  the  Kutta 
condition  at  x  =  1  via  the  relation 

Urn  (1  -x),/J[a*(x) -a'(x)]=0 .  (22) 

A  three-point  Lagrange  extrapolation  formula  is  used 
to  evaluate  the  limit  in  terms  of  A  j.  The  last  equation 
is  obtained  by  imposing  (2)  by  means  of  a  three-point 
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FIG.  2.  Computed  profiles  of  a  sail  with  O'  =ir/4  of  X  =  0.5, 

0.  8,  and  1.0.  The  corresponding  values  of  0  are,  respectively: 
6,  2.  8.  and  1. 6.  The  unit  length  in  which  both  x  and  y  are 
measured  is  c. 

Lagrange  formula. 

For  given  values  of  x  and  a  this  system  of  equations 
is  solved  by  Newton’s  method.  The  values  of  0  andC* 
can  then  be  computed  by  quadrature  from  (16)  and  (17). 

IV.  DISCUSSION  OF  THE  RESULTS 

We  have  used  the  method  described  in  Sec.  in  to 
compute  solutions  for  various  values  of  x  and  a .  Many 
different  families  of  solutions  exist.  They  have  been 
studied  for  a  small  by  Thwaites.1  These  families  can 
be  computed  for  arbitrary  a  and  x  in  the  following  way. 
For  a  given  value  of  x  we  compute  the  solution  for  a 
small  value  of  a  by  using  the  uniform  stream  (i.e.,  the 
solution  corresponding  to  a  =0)  as  the  initial  guess. 
Provided  that  x  is  not  close  to  any  of  the  critical  values 
described  by  Thwaites,1  the  numerical  scheme  conver¬ 
ges  rapidly  to  his  solution  for  the  given  value  of  x. 

Once  a  given  solution  has  been  obtained,  however,  a 
type  of  “boot-strapf’  technique  is  employed,  that  is,  a 
converged  solution  for  one  value  of  X  is  used  as  the  ini¬ 
tial  guess  for  a  solution  with  x  altered  by  a  few  percent. 
Similarly,  tor  fixed  x,  the  angle  a  could  be  varied 


FIG.  3.  Values  of  the  Weber  number  X  aa  a  function  of  0  for 
n  =*/♦.  The  dashed  curve  corresponds  to  Vanden-Broeck  and 
Keller’s  asymptotic  formula. 


FIG.  4.  Values  of  the  lift  coefficient  CL  as  a  function  of  0  for 
a  =  tt/ 4 ,  The  dashed  curve  corresponds  to  Vanden-Broeck  and 
Keller’s  asymptotic  formula. 

using  a  neighboring  solution  as  a  starting  guess. 

Among  all  these  possible  solutions,  only  one  is  “  rea¬ 
sonable”  in  the  sense  that  it  has  a  simple  unimodal 
shape  cambered  to  lee  of  the  wind.  Tuck  and  Hasel- 
grove’s6  stability  analysis  shows  that  this  family  of 
solutions  is  the  only  stable  one  for  a  small.  There¬ 
fore,  we  extend  only  this  family  to  the  nonlinear  case. 

Typical  profiles  are  shown  in  Figs.  1  and  2.  The  val¬ 
ues  of  0  andCL  for  a  =  u/ 4  are  shown  as  functions  of  x 
in  Figs.  3  and  4.  The  broken  lines  correspond  to  the 
asymptotic  formulas  (14)  and  (20)  given  by  Vanden- 
Broeck  and  Keller.5  The  asymptotic  formula  for  0 
agrees  with  the  numerical  results  within  fifteen  percent 
for  x  <  0.3. 

In  Fig.  5  we  present  the  values  of  x  vs  0  for  a  =n/100 
and  a  =w/10.  The  broken  line  corresponds  to  Thwaites’ 
solution.  It  shows  that  his  solution  is  a  good  approxi¬ 
mation  for  a  small,  but  it  is  not  uniform  as  0  tends  to 
zero.  The  discrepancy  between  his  solution  and  the 
exact  numerical  solution  increases  rapidly  as  0  gets 
close  to  zero.  In  the  particular  case  0  =  1,  the  values 
of  C L  were  computed  for  various  values  of  a.  The  nu¬ 
merical  values  of  C L  were  found  to  agree  with  Thwaites’ 
solution  within  5%,  for  a  s  n/18,  20%  for  a<  ir/6,  and 
30%  for  ««  n/4. 


x 


FIG.  5.  Values  of  the  Weber  Number  X  as  a  function  of  0  for 
°  =1/100  (curve  a)  and  a  =  a/1  o  (curve  b).  The  dashed  curve 
corresponds  to  Thwaites'  calculations. 
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Figures  3  and  5  and  similar  results  obtained  for  dif¬ 
ferent  values  of  or  indicate  that  for  each  value  of  a 
there  exists  a  maximum  value  \0(a)  of  the  Weber  num¬ 
ber  above  which  there  is  no  solution  of  the  kind  con¬ 
sidered  here.  For  a  given  value  of  a  there  are  two 
solutions  for  each  value  of  the  Weber  number  in  an 
interval  below  the  maximum  value.  This  nonuniqueness 
is  likely  to  be  of  mathematical  interest  only.  The  rea¬ 
son  for  this  is  that  the  branch  on  the  left  of  the  maxi¬ 
mum  is  probably  unstable.  After  all,  it  corresponds  to 
an  increase  of  the  tension  arising  from  an  increase  in 
the  slack.  This  suggests  that  the  physically  realizable 
solutions  corresponds  to  X  »  X0 (a). 

As  a  increases,  the  slope  of  the  profile  increases  at 
the  leading  edge.  This  was  found  to  limit  the  accuracy 
of  the  numerical  scheme.  Accurate  solutions  for  a 
s  tt/3  and  x  ~X0(o)  could  not  be  computed  with  N  ^  35. 

Finally,  it  is  worthwhile  mentioning  that  the  present 
model  is  not  physical  for  a  large  since  separation  is 


then  likely  to  occur.  A  more  realistic  model  for  a 
large  was  proposed  by  Dugan.4 
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